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Structural quantum criticality and superconductivity 
in iron-based superconductor Ba(Fei_jj:COj^:)2As2 
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We investigated the elastic properties of the iron-based superconductor Ba(Fei_;tC0;t)2As2 with eight Co 
concentrations. The elastic constant C^^ shows large elastic softening associated with the structural phase 
transition. The C^e was analyzed base on localized and itinerant pictures of Fe-3d electrons, which shows 
the strong electron-lattice coupling and a possible mass enhancement in this system. The results resemble 
those of unconventional superconductors, where the properties of the system are governed by the quantum 
fluctuations associated with the zero-temperature critical point of the long-range order; namely, the quan- 
tum critical point (QCP). In this system, the inverse of C^e behaves just like the magnetic susceptibility in 
the magnetic QCP systems. While the QCPs of these existing superconductors are all ascribed to antifer- 
romagnetism, our systematic studies on the canonical iron-based superconductor Ba(Fei_;tC0;t)2As2 have 
revealed that there is a signature of "structural quantum criticality" in this material, which is so far without 
precedent. The elastic constant anomaly is suggested to concern with the emergence of superconductiv- 
ity. These results highlight the strong electron-lattice coupling and eff'ect of the band in this system, thus 
challenging the prevailing scenarios that focus on the role of the iron 3d-orbitals. 

KEYWORDS: iron-based superconductor, elastic constant, quantum criticality, electron-lattice coupling 



1. Introduction 

Iron-based superconductors have been studied actively 
worldwide as new high-temperature superconductors. Their 
superconducting transition temperature is as high as that 
of cuprates, and superconductivity emerges in the system 
with the ferromagnetic Fe element. These characteristics 
are evoking global interest in the mechanism of supercon- 
ductivity.^^ Among the many iron-based superconductors, 
Ba(Fei_;cCo;c)2As2 is suitable for basic research, because it 
provides high-quality large single crystals. The crystal struc- 
ture of the parent compound, BaFe2As2, is shown in Fig. 
l(a).^^ It has a tetragonal crystal structure at room tempera- 
ture. This becomes an orthorhombic structure at the structural 
transition temperature Ts = 140 K accompanying the appear- 
ance of the long-range antiferromagnetic order^^ as shown in 
Fig. 1(b), which is viewed from the c-axis.^^ Ts decreases 
when Co is substituted for Fe, and superconductivity appears 
at X = 0.03."^'^^ Ts falls to zero at x = 0.07, where the super- 
conducting temperature T^c becomes highest, reaching 25 K. 
The magnetic ordering temperature T^ coincides with Ts for 
X = 0, and for samples in which x is less than 0.03, T^ is 
always lower than Ts. Such a phase diagram showing the co- 
existence of magnetic order and superconductivity has been 
observed in the rare-earth compound CePd2Si2 and the ura- 
nium compound UGe2, where superconductivity appears near 
the quantum critical point (QCP).^'^^ 

Intensive studies on QCP have led to a new paradigm in 
physics in which magnetic fluctuation mediates superconduc- 
tivity. In the case of iron-based superconductors, the mech- 
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anism of superconductivity has so far been discussed based 
on the origin of the adjacent magnetic and structural phases, 
whether magnetic or orbital.^' Mazin et al. and Kuroki et 
al. predicted that a fully gapped sign-reversing s-wave state 
(s+ state) would be realized by magnetic instability.^^' On 
the other hand, Kontani and Onari, and Yanagi et al. pro- 
posed a conventional s-wave state without sign reversal (s++ 
state) mediated by orbital fluctuation.^^' Investigations on 
the neighboring order and its fluctuation are important to clar- 
ify the superconducting properties of iron-based supercon- 
ductors. 

Our experimental tool, ultrasonic measurements, provides 
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Fig. 1. (Color online) (a) Crystal structure of BaFe2As2, which belongs to 
the base-centered tetragonal crystal class h/mmm, and (b) crystal structure 
and magnetic order below Ts with the crystal symmetry of Fmrnm- 
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US with information on changes in the symmetry of the lat- 
tice system. Strains introduced into the crystal in the elas- 
tic constant measurements deform the crystal locally and 
break the crystal symmetry (symmetry breaking field). When 
the system encounters either a charge, or orbital (electric 
quadrupole) order, the elastic stiff'ness that possesses the same 
symmetry as, and thus couples with, this order tends to exhibit 
anomaly (softening in most cases). Accordingly, by examin- 
ing the temperature (r)-dependence of the anisotropic elastic 
stiff'ness Cij, we can obtain fundamental information on the 
symmetry of the order. 

The T-dependence of the elastic stiff'ness of 
Ba(Fei-xCox)2As2 has been reported for platy samples 
of BaFe2As2 and Ba(Feo.8Coo.2)2As2 using the resonant ul- 
trasonic spectroscopy (RUS) method^"^^ and for a bulk sample 
of overdoped Ba(Feo.9Coo.i)2As2 using pulsed ultrasonic 
spectroscopy with a phase comparison method. The former 
study reported a large elastic softening toward Ts in Cee, 
which was present at the percentage ratios of 95% and 18% 
in BaFe2As2 and Ba(Feo.8Coo.2)2As2, respectively. The latter 
study reported that Cn, C33, C44, and j (Cu - Cu) gradu- 
ally increased as T decreased, and showed no remarkable 
anomaly at T^c in contrast with the Cee anomaly of 21%. 

In the present study, we focused our attention on the 
huge anomaly in Cee and measured the T-dependence of 
Ba(Fei_;cCO;c)2As2 with eight Co concentrations from x = 
to X = 0.245 to investigate the neighboring phase and its role 
in the superconducting properties. 

2. Experimental 

2. 1 Ultrasonic measurements 

We measured the elastic constants for Ba(Fei_;cC0;c)2As2 
using an ultrasonic pulse-echo phase comparison method 
as a function of temperature from 5 K to 300 K using a cryo- 
stat mounted on a Giff'ord-McMahon (GM) cryocooler. Ul- 
trasound was emitted and detected by LiNbOs transducers. 
X-cut plate of LiNbOs at 41° with a thickness of 100 yum was 
used for the transverse waves. The fundamental resonance fre- 
quency was 19 MHz. We usually use their higher harmonics 
of 67 MHz for the experiments. The transducers were glued 
onto a pair of parallel planes of single crystals using an elastic 
polymer Thiokol. 

Elastic stiff'ness was obtained by C = pv^, where p is 
the density and v is either the longitudinal or transverse 
sound velocity. The corresponding sound velocities can be 
obtained by choosing the propagation and displacement direc- 
tions. Tetragonal crystal symmetry has six independent Q/s; 
namely, Cu, C33, C12, C13, C44, and C^e- The propagation and 
displacement directions of the sound velocity was [100] and 
[010] for Cee- The sound velocity was obtained by the time 
interval of the echo train and the sample length, whose accu- 
racy is within a few percent due to the usage of large crystals. 
The value of p was obtained using the data of x = and 0. 1 ."^^ 
By assuming Vegard's law, the lattice constant of the a(b) and 
c axes 8iTea = b = 0.39636 -h 3.8981 10"^ x (nm) and c = 
1.3022 - 0.0421JC (nm), respectively, for Ba(Fei_;,Co;,)2As2. 
The value of p is 6.48 x 10^ kg-m"^ and 6.55 10^ kg-m"^ for x 
= and 0.245, respectively. The XYZ coordinate was defined 
by the unit cell of the h/mmm crystal structure. To prevent 
damage to the sample due to rapid changes in temperature. 
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Fig. 2. (Color online) Temperature-dependence of the elastic stiffness C66 
of Ba(Fei_;cCO;c)2As2 with various values. 

the rate of change in temperature was carefully controlled so 
as to be 10 K/h near T^?^ 

High-quality large single crystals of Ba(Fei-xCox)2As2 
used in this work were grown by the self-flux method. Sam- 
ples with eight Co concentrations x = 0, 0.037, 0.060, 0.084, 
0.098, 0.116, 0.161, and 0.245 were prepared. The Co con- 
centration in the grown crystals was determined by energy- 
dispersive X-ray spectroscopy (EDS) measurement. The Co 
content of the samples was about 75% of the prepared one. 
The samples were cut into a rectangular shape, after determin- 
ing their axis by X-ray Laue photography. The typical size of 
a sample was 3 mm x 3 mm in the tetragonal ab (XY) cleav- 
age plane, and 2 mm in thickness on the c(Z)-axis. 

2.2 Temperature dependence ofC^e 

Figure 2 shows the T-dependence of Cee for samples with 
x = 0, 0.037, 0.060, 0.084, 0.098, 0.116, 0.161, and 0.245. 
Cee significantly decreases as T decreases. While a decrease 
in elastic stiff'ness is a common feature as a precursor of the 
structural phase transition, the amount of softening, as much 
as 80% for Ba(Feo.963Coo.o37)2As2, is unprecedentedly large. 
The softening in C66 corresponds to the symmetry change 
from tetragonal to orthorhombic. This is consistent with the 
structural analysis of this material, where the space group is 
h/mmm and Fmmm for the high- and low-temperature phases, 
respectively.^^ 

The decrease in C66 with decreasing T is prominent for x 
= 0, 0.037, and 0.060, which follow the disappearance of the 
signal at the structural transition temperature Ts = 141 K (x 
= 0), 84.7 K(x = 0.037), and 30 K (jc = 0.060). The data 
below Ts were not plotted for x = and 0.037, because the 
sound echo signal disappeared at a certain temperature range 
below Ts, which may be ascribed to strong scattering of the 
sound by the orthorhombic domain boundaries. The softening 
of C66 is less prominent as x increases, which is consistent 
with the disappearance of the structural phase transition. For 
the superconducting samples, 0.060 < x < 0.116, anomalies 
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Fig. 3. (Color online) Temperature-dependence of the inverse of S^^^cr- 

are observed at T,^ = 24.0 K (jc = 0.060), 20.7 K (jc = 0.084), 
16.5 K (jc = 0.098), and 10.5 K (jc = 0.116). 

The underdoped samples show a peak at T^c and a step-Hke 
anomaly below T^^, as shown in Fig. 8. On the other hand, 
Cee increases below T^^ in the overdoped region. This differ- 
ence will be discussed later. Moreover, the superconducting 
transition disappears when x > 0.161. C^e for the samples of 
X = 0.161 and 0.245 shows gradual increases as T decreases, 
which a typical behavior for any material, reflecting phonon 
anharmonicity. 

3. Discussion 

3. 1 Analysis based on localized picture 

The observed large elastic anomaly in C66 is a precursor 
of a structural phase transition. Here, we will discuss the ori- 
gin of the softening in C66 from the viewpoint of experimen- 
talists to provide the analysis method. ^'^^ The 3d orbitals are 
split into Eg doublet and T2g triplet by the electric crystalline 
field (CEF) in a cubic symmetry. The Eg is spilt into two sin- 
glets, and the T2g is split into one singlet and one doublet by 
tetragonal CEF of iron-based materials. The remaining dou- 
blet can be lifted by the elastic strain, and may cause the elas- 
tic anomaly. 

We will consider the coupling between the strain s and an 
order parameter O sls H = -AOrSr. The equivalence of the X 
and Y axes in tetragonal symmetry leads to the degeneration 
of the dzx and dyz orbitals. This degeneracy is lifted by sxy 
or Exx - ^YY, and brings about anomalies in the correspond- 
ing Cee and ^ (Cn - Cu). The strain sxy can couple with the 
orbital (quadrupole) Oxy- We can calculate the strain suscep- 
tibility based on the localized picture of d electron. 

C66 = C66,0 - NA^- ^-r- = C66,0— 7^ (1) 

where A, /, and N are the coupling constant, the intersite in- 
teraction, and the number of atoms per unit volume. Here, we 

adopted the form ofx^^ = ^^^^j^ ^ , for localized d electrons. 



Here, we introduce the elastic compliance Sij, which is a 
component of the inverse Cij matrix. In Eq. (1), the transi- 
tion undergoes at Tq, where the lattice shows an instability. 
Elastic compliance represents the gstructuralh susceptibility 
of elastic systems, and corresponds to the magnetic suscep- 
tibility X in magnetic systems. The experimentally observed 
^66 (= l/Qe) can be decomposed into the sum of the anoma- 
lous contributions that exhibit critical behavior S e6,cr and the 
normal contribution (background) 5 66,o- 

See = Sl^-\- S 66,cr = S 66,0 |l + YZJ^ 

where Ejj = Tc - 0. Ejj stands for the Jahn-Teller energy, an 
energy scale that corresponds to the strength of the electron- 
lattice coupling. Note that this formula has the same form as 
the Curie- Weiss susceptibility of ferromagnetic materials. 

For the analysis, we employed the data on 
Ba(Feo.755Coo.245)2As2 for 5'66,o, and subtract it from 
the other data. Figure 3 shows the inverse of S e6,cr as a func- 
tion of temperature. It can clearly be seen that 1/Se6,cr in the 
underdoped region (x < 0.070) exhibits linear T-dependence. 
This indicates that 5'66,cr obeys Eq. (2). 

In the underdoped samples, Tq in Eq. (2) and Ts obtained 
by C66 are 141 K and 136 K for x = 0, 84 K and 81.4 K for x 
= 0.037, and 30 K and 28.2 K for x = 0.060, respectively. The 
closeness of the values in each case suggests the occurrence 
of ferro-order. This strongly suggests the crystal deformation 
of Sxy below Ts. 

For the underdoped samples, E]j is approximately equal 
to 50 K. This value is comparable to A = 20 K in Ref. 
[15], which is smaller than our value. The diff'erence between 
the two experiments is originated from how to choose the 
background Cq. While we analyzed the data without a fit- 
ting parameter for the background by adopting the data of 
Ba(Feo.755Coo.245)2As2 for Co, Goto et al. carried out the fit- 
ting including the background as adjustable parameters. We 
performed the analysis by using the same Co, which can re- 
veal the systematic behavior though the whole Co concen- 
tration. The values of the electron-lattice coupling constant 

^ -yjiiOxY)^) was calculated from Ejj, and is 0.25 eV/Fe for 
X = 0.037. 

3.2 Analysis based on band picture 

When X exceeds 0.07, 5'66,cr deviates from T-linear behav- 
ior below a certain temperature T*; i.e., 40 K for x = 0.084, 
130 K for X = 0.098, and 190 K for jc = 0.116. The T- 
dependence of C66 in the overdoped region can be explained 
by a band picture instead of Eq. (2). Large elastic anomalies 
compared with iron-based materials have been reported in the 
A15 superconductor V3 Si and the Laves-phase superconduc- 
tor CeRu2 so far.^^' These anomalies have been ascribed to 
the large density of states at the Fermi energy. The 3d-orbitals 
form bands in an iron-based superconductor. The bands lo- 
cated above the Fermi energy at the -point form hole Fermi 
surfaces and electron pockets at M-points of the zone bound- 
aj^y 13) Yom M-points (which are identical to X-points in the 
square lattice formed by Fe atoms) do not becomes equivalent 
under exy-^^^ If a large density of states exists at the M-point, 
an elastic anomaly may be caused. This provides a natural so- 
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Fig. 4. (Color online) Theoretical fitting of C66 by the band Jahn-Teller 
effect under the assumption of Lorentzian density of states. Adjustable pa- 
rameters are the bandwidth W and the intersite interaction /. 
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lution to the question of why the anomaly only emerged in 
C66.^^^ The formula for the elastic constant based on this con- 
sideration is as follows:^^'^^^ 



C = Co - (dMi - dui) 



(3) 



xl=i^i:fk{\-fk) = SdEN{E)f{E) 



keT 

N(E) = ^ 



1-1 



where / is the Fermi-Dirac function, d (= Jmi = -dMi) is the 
electron-lattice coupling constant, Nq is the density of states 
at Fermi energy E^, and / is the intersite interaction. The re- 
sults of fitting by adopting the Lorentzian density of states are 
shown in Fig. 4. 

In this analysis, the main adjustable parameter is the band- 
width W. We fixed Ey^/W = 0.4 to obtain the best fitting. The 
value of / is displayed in the figure, depending on the value 
of X. We added a small constant ranging from 0.8 GPa to 
2.45 GPa to the experimental value of Co for the fitting. The 
value of W was 100, 180, 250, and 450 K for the 8.4%, 9.8%, 
11.6%, and 16.1% samples, respectively. From the adjustable 
parameter Nod^, the value of d was respectively evaluated as 
0.22, 0.25, 0.28, and 0.28 eV/Fe for the four samples, under 
the assumption that A^o is approximately equal to l/W. The 
values of A^o calculated from Nq = l/W are respectively 110, 
64, 46, and 26 states/eV for the four samples. 

The band eff'ect discussed here is called the electronic redis- 
tribution mechanism,^^^ and is applicable to systems having 
low resistivity. It is supported experimentally. The resistivity 
of this system becomes lower in the overdoped region com- 
pared with the underdoped region. ^^'^"^^ Band nesting along 
the [7i/a,7i/a,0] direction is a characteristic feature of iron- 
based superconductors. We have no tool for analyzing elas- 
tic data based on band nesting. However, we infer that it has 



a similar eff'ect to that discussed above, and may cause C66 
softening. The obtained parameters are summarized in Fig. 5 

3.3 Quantum criticality 
3. 3. 1 Phase diagram 

Figure 6 summarizes the phase diagram of the 
Ba(Fei_;cCo;c)2As2. The Co concentration-dependence 
of Ts and Tsc are highly consistent with the values reported 
elsewhere. ^^'^^^ We found two characteristic temperatures of 
r* and Tmax. A possible explanation of is the crossover 
from the non-Fermi liquid to the Fermi liquid region. The 
boundary from non-Fermi T-dependence to Fermi liquid 

-dependence observed in the resistivity measurements is 
not clear for Ba(Fei_;cCO;c)2As2.^^^ However, the behavior 
observed in BaFe2(Asi_;cPx)2 is similar to the case in the 
present study including the x-dependence of the crossover 
region. 

We found r^ax, which corresponds to the temperature at 
which 5'66,cr takes the maximum value (l/5'66,cr takes the min- 
imum value), as shown in Figs. 3 and 7. For highly corre- 
lated electron systems such as CeCu2Si2, UPd2Al3, and UPts, 
a similar maximum was reported in the magnetic susceptibil- 
ity;^, and is interpreted as a Kondo temperature, signaling the 
coherent motion of f-electrons.^^^ From this analogy, r^ax in 
this system may suggest the onset of new coherent states. 

It would be notable that T* and Tmax go to zero at the QCP 
concentration of xc = 0.07, indicating the existence of struc- 
tural QCP at this concentration. In addition, the analysis based 
on the band picture suggests that the bandwidth W also goes 
to zero at the same xc, which suggest a possible mass en- 
hancement toward QCP in this system. Note that this point is 
located at the center of the superconducting dome. A similar 
phenomenon was reported for BaFe2(Asi_;cPx)2.^^^ The ob- 
tained phase diagram and various phenomena near the QCP 
resemble those of the well-known rare-earth compounds and 
uranium compounds. This coincidence strongly suggests the 
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Fig. 6. (Color online) Phase diagram of Ba(Fei_;tCO;t)2As2. Filled circles: 
Ts, and filled triangles: T^c from other studies. ^^'^^^ Open symbols in this 
study. The curves are the guides for eyes. 



intimate relationship between superconductivity and QCP in 
this system, similarly to heavy fermion systems where mag- 
netic QCP is believed to be responsible for the emergence of 
unconventional superconductivity.^^ The essential difference 
in this case is that the quantum criticality is associated with 
structural fluctuations, unlike the magnetic fluctuations in the 
previous cases. From this reason, we would like to name the 
structural quantum criticality. 

3.3.2 Correlation between elastic anomaly and supercon- 
ductivity 

We now discuss the relationship between the elastic 
anomaly and superconductivity. As seen in the inset of Fig. 
7, the amount of l/5'66,cr is proportional iox-xc, where xq is 
the QCP concentration of Co; xc = 0.07 for this system. Such 
behavior is well known in;^ of the magnetic QCP. It is surpris- 
ing that such well-known behavior holds in this system with 
the respect of S 66 instead of x for the magnetic system. As 
shown in the same figure, T^c decreases with increasing x-xq- 
So we can recognize an apparent correlation between T^c and 
l/^66,cr that Tsc is a function of l/5'66,cr The explanation for 
this interesting fact is speculated to be as follows. As shown 
in Fig. 8, the underdoped sample exhibits a small anomaly 
at Tsc, while a large upturn at T^c is seen in the overdoped 
samples. Once the system enters the orthorhombic phase from 
the tetragonal phase, structural fluctuations are suppressed in 
the ordered phase. In the overdoped samples, however, struc- 
tural fluctuations still survive even at Tsc- The amount of the 
anomaly at T^^ correlates with the peak height of 5'66,cr, which 
is the measure of structural fluctuation. The large anomaly at 
Tsc for the overdoped samples suggests strong coupling be- 
tween the structural fluctuations and superconductivity. These 
facts suggest that the origin of S e6,cr is deeply related to the 
emergence of superconductivity. 

4. Summary and concluding remarks 

The microscopic origin of C66 softening has been discussed 
based on various mechanisms. Large elastic anomalies are of- 
ten observed in the 4f electron system with a bilinear coupling 
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Fig. 7. (Color online) Temperature dependence of Se6,cr for the overdoped 
samples. Inset is T^c from this work and Ref. [25], and the inverse of the 
peak value of Se6,cr as a function of the distance from the QCP concentra- 
tion X - XQ. The curve of T^c is the guide for eyes. 



between quadrupole operator O and elastic strain s with the 
form of Os. First, the effect of spin-nematic order has been 
considered, where the elastic anomaly above is caused by 
the fluctuation of a pair of magnons at two sublattices.^"^'^^^ 
Spin-nematic order parameter can couple to the strain with 
a bilinear form. In this case, the nematic order deforms the 
crystal through the spin-orbit interaction. 

On the other hand, theories based on orbital fluctuation 
have been proposed. Yanagi et al. considered Ferro-orbital 
order accompanying large elastic softening. Kontani et al. 
have discussed a two-orbiton process, which is initiated by 
antiferroquadrupolar fluctuation caused by interband nesting 
and consisting of a coupling between optical phonon and or- 
biton at the zone-boundary. They showed that two-orbiton 
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Fig. 8. (Color online) Temperature dependence of C^^ for the overdoped 
Ba(Feo.94oCoo.o6o)2As2 near T^c in expanded scales 
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process brings about very large C^e softening. 

As regards the values of the electron-lattice coupling con- 
stants A, the spin-nematic theory and the orbital theory give 
17 meV/Fe and 0.2 eV/Fe, respectively. Our experiment 
provides the value A ~ 0.22-0.25 eV/Fe for the underdoped 
samples, and d ~ 0.2-0.28 eV/Fe for the overdoped samples. 
These values of the coupling constant are considered to be 
very large, and are consistent with the orbital-based theory 
both from the localized and band pictures. Generally speak- 
ing, a coupling between the magnetic moment and the stain 
is mediated by a spin-orbit interaction. It is considered that 
such strong electron-lattice coupling is not caused by weak 
spin-orbital coupling in transition metals, but by the orbital 
nature of this system. We suppose that the strong electron- 
lattice coupling is a characteristic feature of multi-band sys- 
tems, possessing orbital degrees of freedom. Our results are 
considered to support the orbital-based theory. 

In this article, structural quantum critical behavior was re- 
ported for Ba(Fei_;cCO;c)2As2. The QCP behavior has been 
also reported by resistivity measurements and NMR of 
Ba(Fei_;,CO;,)2As2.^^'^'^^ The QCP behavior observed in NMR 
and elastic measurements would be expected to be ascribed 
to the same origin. Elastic constant is no sensitive probe for 
magnetism, but sensitive to orbital (quadrupole), contrary to 
NMR. Our measurement is the observation from the side of 
orbital (quadrupole). 

Further discussions will be necessary on the role of the spin 
and orbital in the mechanism of superconductivity. However, 
this system was found to be well characterized by a strong 
electron-lattice coupling. Therefore, irrespective of whether 
the spin or orbital mediate superconductivity, the results of 
the present study suggest that the structural fluctuation must 
be actively incorporated into their role, and that the orbital 
degrees of Fe-3d must also be taken into consideration in or- 
der to understand the full picture of superconductivity in iron- 
based compounds. 
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